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Abstract 

Given a homogeneous Poisson point process in M. d , Haggstrom and 
Meester |2] asked whether it is possible to place spheres (of differing 
radii) centred at the points, in a translation-invariant way, so that 
the spheres do not overlap but there is an unbounded component of 
touching spheres. We prove that the answer is yes in sufficiently high 
dimension. 

1 Introduction 

A sphere process is a simple point process A on l d x [0, oo). The support 
of A is the random set [A] = {(x, r) eR d x [0,oo) : A({(x,r)}) = 1}. If 
(x,r) £ [A], we say that there is a sphere of radius r £ [0, oo) at x £ M. d . 
The centre process A is the point process on M. d given by the projection 
A(-) = A(- x [0, oo)). We say that A is a Poisson sphere process if A is a 
homogeneous Poisson process. 

A hard sphere process is one in which the interiors of the spheres do 
not overlap; that is, almost surely 

p{x, y) > r + s for any distinct (x, r), (y, s) £ [A], 
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where p denotes Euclidean distance. (Note in particular that, since all the 
radii are non-negative, in a hard sphere process A no sphere may contain any 
points of A other than its own centre.) For z G M. d the shifted sphere process 
A + z is defined by (A + z)(A) = A(A — z), where A — z := {(x — z,r) : 
(x,r) G A}. The sphere process A is (translation-)invariant if A + z and A 
are equal in law for all z G M. d . 

Let G(A) C M. d be the random set covered by all the spheres: 

G(A) = {y G R d : p(y,x) < r for some (x,r) G [A]}. 

The connected components of G(A) are called clusters. We say that A 
percolates if there is an unbounded cluster. Our main result is the following. 

Theorem 1 For all d > 45, there exists an invariant Poisson hard sphere 
process A which percolates almost surely. 

Haggstrom and Meester j2] studied several invariant continuum perco- 
lation processes, and proved that the "dynamic lily-pond model" does not 
percolate. This is a Poisson hard sphere process in which spheres grow from 
all Poisson points at the same rate, and whenever two spheres touch, they 
both stop growing. Haggstrom and Meester asked whether there exists an 
invariant Poisson hard sphere process which percolates. 

In dimension d = 1 it is easy to see that any Poisson hard sphere process 
almost surely does not percolate. For consider 4 consecutive points x\ < X2 < 
x 3 < X4 of the Poisson process which satisfy X3—X2 > (22 — xi) + (x^— 23); then 
the spheres centred at X2 and X3 cannot touch. Since such a configuration of 4 
points appears infinitely often in a Poisson process, no percolation is possible. 
It is unknown whether there exists a percolating Poisson hard sphere process 
in dimensions 2 < d < 44 (even without the requirement of invariance). 

Our explicit bound of 45 in Theorem [T] could probably be reduced with 
some effort. Our construction certainly cannot be adapted to work in dimen- 
sion less than 3 (and probably not in dimension that low). 

Jonasson [3] showed that in the hyperbolic plane H 2 there exists a hard 
sphere process that percolates when the Poisson process has sufficiently low 
intensity, but was unable to determine what happens for high intensity. (In 
M. d , scaling shows that the intensity of the process is immaterial). 

Note that our definition of a hard sphere process allows spheres of radius 
zero. Our proof of Theorem ^ can in fact be adapted to prove the existence 
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of an invariant percolating Poisson hard sphere process in which the spheres 
all have positive radii. We explain this in a remark at the end of Section 0] 

Our proof of Theorem^is in two parts. First we construct a non-invariant 
hard sphere process T which percolates. Then we convert this to an invariant 
process A by "stationarizing" - applying a uniform random translation in a 
large ball, and taking a limit. The non- invariant construction of T proceeds as 
follows. Starting from a Poisson process, we attempt to grow an unbounded 
cluster iteratively. At each step we try to choose a radius for a new sphere 
(centred at some Poisson point) so that it touches the cluster constructed so 
far. This will be possible only if the proposed new sphere contains no other 
Poisson points, and we need to ensure this happens sufficiently often that 
the cluster can continue to grow. We do this by comparison with a certain 
two-dimensional percolation process; in sufficiently high dimension we can 
arrange that the probability of success at each step exceeds the relevant 
critical probability. 

To ensure that the stationarized version A also percolates, the unbound- 
ed clusters of T should occupy a positive fraction of space - see below 
for a precise statement. We will achieve this by repeating the essentially 
two-dimensional construction described above in infinitely many "layers" 
throughout R d . 

Let L denote Lebesgue measure on M. d , and denote the Euclidean ball 
B(x, r) = B d (x, r) := {y e M. d : p(x, y) < r}, and the origin = (0, . . . , 0). 
Define the lower density of a set A C M. d to be 



Using the construction sketched above, we shall prove the following result, 
from which Theorem Q will be deduced. 

Theorem 2 For d > 45 there exists a (not necessarily invariant) Poisson 
hard sphere process T which percolates almost surely. Furthermore, T can be 
chosen so as to have the following additional properties. 

(i) The union I of all unbounded clusters satisfies lim^oo EV(I^) = 1. 



V(A) = liminf 



C(A(lB(0,r)) 



r— >oo 




For s > define the s-neighbourhood of a set A C M. d to be 



= |J B(x,s). 
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Figure 1: Part of the lattice H* 

(ii) There exists a constant K(d) < oo such that T(M. d x (K, oo)) = 
almost surely; that is, there are no spheres of radius greater than K. 

The article is organized as follows. In Section |2 we give the construction 
that gives rise to the non-invariant process in Theorem |21 In Section |3] we 
prove that the resulting process does indeed have the properties stated in 
Theorem |21 In Section 0] we deduce Theorem ^ 

2 Construction 

In this section we describe the construction of the hard sphere process in 
Theorem |21 We will prove that it has the required properties in the next 
section. The construction is given in terms of parameters fi = 0.75, 5 = 0.1, 
e = 0.01 and C, where C = C(d) is a (large) constant to be chosen later; 
the values for these parameters are chosen so that the construction yields an 
unbounded cluster. Let A = \(d) > be another constant to be chosen later, 
and let II be a homogeneous Poisson process of intensity A, and denote its 
support [II]. We will construct a hard sphere process T such that T = II. 

Let EI be the hexagonal lattice in IR 2 , with edge-length 2. (Thus, the faces 
are regular hexagons of side 2, and there is a vertex at the origin = (0, 0), 
say.) Let H* be the graph formed from EI by adding an extra vertex in the 
middle of each edge; see Figure ^ We will call the vertices of the original 
lattice H site vertices and the extra vertices of H* bond vertices. Also 
fix some arbitrary well-ordering of the vertex set V(H*) of M*. 

Fix the dimension d > 3. For any vertex v G V^H*), define the cell 
centred at v to be the set 

W(v) := B 2 (v, e) x B d _ 2 {0, C) C R d . 
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Figure 2: An example of the construction in dimension d = 3. Shown are 
the cells centred at the origin and at the three neighbouring bond-vertices 
u,v,w, together with two touching balls B(x ,r ) and B(x u ,r u ). 

See Figure 121 

First, here is a brief description of the construction. We will attempt to 
place spheres with their centres in distinct cells in such a way that spheres 
at adjacent cells touch; see Figure El All the spheres will have radii in 
\p — 5, \i + 5] . Note that any two non-adjacent vertices of HI* are at distance 
at least ^/3. Since 2(/x + 5 + e) = 1.72 < y/3, it will be impossible for spheres 
centred in non-adjacent cells to touch or overlap. 

We will explore the lattice iteratively starting from the origin, attempting 
to construct an unbounded cluster. Each step of this exploration will have 
two parts. Firstly, for some vertex w G V(H*) we try to find a point y w G 
[II] Pi W(w) such that it is possible to place a sphere of some radius r w G 
[/i — S, n + 6] at y w that touches one of the existing spheres and does not 
overlap any of the existing spheres. Secondly, we check to see whether or not 
there are any other points of the Poisson process within B(y w ,r w ). If there 
is none, our construction succeeds and we let x w = y w . If either part fails, 
we let x w = A. 

At each step of the construction, each vertex v G V(H*) will be either 
good (meaning we succeeded in constructing a sphere B(x v ,r v ) at v), bad 



(meaning the construction failed), or unexplored. 

Using HI* (instead of H, say) will enable us to explore in such a way that 
every time we explore a new vertex it is adjacent to exactly one good vertex 
and no bad vertices. This will simplify our arguments. Other choices of 
lattice are possible, such as the square lattice in place of HI - see the remark 
at the end of Section El 

Here is a formal description of the construction. As remarked earlier, we 
will try to construct a cluster in each of a set of two-dimensional "layers" . 

First Layer 

Start with all vertices in V(H*) unexplored. We perform a sequence of steps 
0, 1, 2, . . .. One new vertex (or sometimes two) will be explored at each step. 

Step 

Let wq = 0; we will start our exploration at the site vertex wq. Temporarily 
denote its cell W = W(wq). Consider two cases: 

Case 1: If W n [II] = 0, then let x = y = A. 

Case 2: If W D [II] 7^ 0, then pick yo uniformly at random from the set 
W n [II] (conditional on W [II]). Take r = fi. Now take 

/ y Q if B(y ,r )n[U] = {y }; 
Xo \ A if B(y ,r )n[U]D{y }. 

Declare wq = to be bad if xq = A, and good otherwise. 

Step n (n > 1) 

Step n consists of parts (a)-(c). 

(a) We first choose an unexplored vertex w n to explore, according to the 
following rules. 

(i) If there is an unexplored site vertex w adjacent to some good bond 
vertex and two unexplored bond vertices, choose w n to be the first 
such w in the ordering on V(M*). 
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(ii) If there is no w as in (i), but there is an unexplored bond vertex 
w adjacent to a good site vertex and an unexplored site vertex, 
choose w n to be the first such w. 

(iii) If neither (i) nor (ii) hold then stop, and proceed to 'Subsequent 
Layers' below. 

(b) Temporarily write w = w n for the vertex chosen in (a). The rules in (a) 
ensure that w is unexplored, and it has exactly one good neighbour, v 
say, while all the other neighbours of w are unexplored. We will try to 
construct x w G W(w), the centre of a sphere tangent to B(x v ,r v ). 

Define 

U(v) = B(x v ,r v + n + 5) \ B(x v ,r v + /i - 8); 

this is the set of possible centres for a sphere of radius G [pi — S, + 5] 
that touches B(x v ,r v ). Let 

s = w(w)nu(v). 

Now consider two cases: 

Case 1: If S n [n] = 0, then let x w = y w = A. 

Case 2: If Sn [II] ^ 0, then pick y w uniformly at random from SO [II] 
(conditional on S H [II]), and let 

(by the definition of V(v) we have r w G [p, — S, \i + 6]). Now 
take 

f y w if B(y w ,r w ) H [II] = {y w }; 
Xw \ A if B(y w ,r w )n[U]D{y w }. 

(c) Declare the vertex w to be bad if x w = A, and good otherwise. 

In addition, if w is a bond vertex and is declared bad, then declare its 
remaining unexplored site vertex neighbour to be bad also. 

Continue either for an infinite sequence of steps, or until we stop in 
(a) (iii) above. 
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Subsequent Layers 

All the points x w constructed above (for good vertices w) lie in the "layer" 
M 2 x l?d_ 2 (0, C), and the radii r w are at most (j, + S. We want to repeat 
the construction in other layers, ensuring that the spheres in different layers 
do not overlap. Therefore let L = 2(C + \i + 8 + 1), and consider the lat- 
tice of points LZ d ~ 2 C M. d ~ 2 . For each z G LZ d ~ 2 , repeat the "first layer" 
construction above, but now in the layer M 2 x B d ^ 2 {z,C) (using the Poisson 
points in cells of the form B^{v, e) x Bd-2{z, C)). This results in independent 
identically distributed clusters in each of the layers. 

Definition of T 

Finally, construct the hard sphere process T by placing a sphere centred at 
x w with radius r w , for each good vertex w, in every layer. Also place a sphere 
of radius zero centred at each remaining Poisson point z G [IX] . 

3 Proof of Construction 

In this section we prove Theorem El using the construction of Section El We 
start by assembling some tools. 

Percolation 

Let EI be the hexagonal lattice in R 2 with side-length 2. Consider Bernoulli 
site percolation with parameter p on HI. That is, each vertex is open with 
probability p and closed otherwise, independently for different vertices. There 
exists a critical probability p c = p^ lte (EI) < 1, with the property that if p > p c 
then there is a.s. a unique infinite connected cluster of open vertices (see [T]). 
It is proved in [7j that 

p c < 0.794. 

Now suppose p > p c and let K C M? be the vertex set of the infinite 
open cluster if contains 0, and let K = otherwise. Suppose d > 3 and let 
(K z ) zeZ d-2 be a family of i.i.d. random sets each with the same law as K. 
Fix any L > 0, and define the random set 

Y = U K z x {Lz} C R d . 
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Lemma 3 For p > j9^ ltc (HI) and any L > 0, the set Y defined above satisfies 
as a — > oo. 

Proof. Let / be the infinite open cluster, and let 6 = P(0 G I) > 0. Define 
a random set as follows. Flip a #-coin independently of J, and let = / 
with probability 9, and K = with probability 1 — 0. By the Harris-FKG 
inequality (see for example £Q) we see that K stochastically dominates K. 
Now let (K z ) zeId d-2 be a family of i.i.d. random sets with the same law as 
K, and let Y = \J zeZ d-2 K z x {Lz} C R d . Then F dominates Y, so D(yW) 
dominates T>{Y^). The random set F^ is invariant in law under isometries 
of HI x Z d_2 , so Fubini's theorem implies 

EV(Y {a} ) > inf P(x G F w ) = inf P(B(x,a) HF ^ 0) -> P(F ^ 0) = 1 
as o — > oo. □ 



Random points 

Let LT be a point process in R d . For r > we call a point x G [II] r-isolated 
if there is no other point of [II] within distance r of x. 

Lemma 4 Let U be a homogeneous Poisson point process of intensity A in 
R d . Let S C R d be a Borel set with CS G (0,oo). Let X be a point cho- 
sen uniformly at random from the set [II] fl S (provided it is non-empty), 
conditional on II. Then 

P (X exists and is r-isolated) > e ~ XCB{Q ' r) - e" A£5 . 

Lemma 5 Under the assumptions of Lemma^Q let Z C R d 6e disjoint from 
S. Then 

P(X exists and is r-isolated | n(Z) = 0) > P(X exists and is r-isolated). 

Proof of Lemma QJ For convenience, let X = A if 11(5') = 0, and take 
A to be not r-isolated. Denote the random sets B = B(X,r), U = B C\S, 
and V = B \ S. Then we have 

P(X is r-isolated | U(S),X) 

= p(n(c/) = i,n(v) = o | n(S), x)i[n(S) > o] 

= (l-^J e-^lpl(5) > 0]. 
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Hence 



P(X is r-isolated) 




E 



e 



—XCB 



□ 



Let U\s denote n restricted to S; that is the point process with support 

[n] n s. 

Proof of Lemma [HI We have 

P(X is r-isolated, U(Z) =0\lL\ s ,X) 
= 1 [U(B(X, R)nS) = l]P (U(B(X, R)\S) = 0, U(Z) = | X) 
> l[U(B(X,R) nS) = 1]P(U(B(X,R) \S) = | X)P(U(Z) = 0) 
= P(X is r-isolated | U\ S ,X)P(U(Z) = 0). 

Taking expectations yields the result. □ 
Volume bound 

We write ujd := £Bd(0, 1) for the volume of the unit ball. Fix d > 3 and 
e, C > 0. Suppose w, w' G M 2 are such that p(w, w') = 1. Define the sets 

W = B 2 (w,e) x B d _ 2 (0,C) and W = B 2 (w', e) x S d _ 2 (0, C). 
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Lemma 6 Fix /i = 0.75, 5 = 0.1 and e = 0.01. Let d> 11. Let W, W be as 
above. There exists C = C'(d) such that if C > C , then for all x G W and 
r G [fi — 5, fi + 5], writing 



S = W'n (B(x,r + fx + 5)\B(x,r + fi-5)), 



we have 



CS > 



d-2 

1.2— - 1 



) 



3 



In order to prove Lemma El we use some further geometric results. 

Lemma 7 Fix d and < R- < R + < oo. There exists C = C'(d, _R_, R + ) 
such that if C > C , for all x G -B(0, C) and R G {R-, R+) we have 



Proof. First suppose that p(0,x) = C (so that x is on the surface of the 



and note that / depends only on the ratio R/C, while for fixed R, the 
function / is increasing and continuous in C, and converges to 1/2 as C — > oo 
(since near x, the ball B(0,C) approaches a half space). Therefore by the 
intermediate value theorem the claimed result holds for the case p(0, x) = C. 

Suppose now that p(0,x) < C. The result is trivial when B(x,R) C 
5(0, C). If not we can replace B{0, C) with 5(0, C), where p(0, x) = C (and 
so C G [C — R + , C)), and appeal to the case already proved. □ 

Lemma 8 Let /i, 5, e and W, W be as in Lemma® There exists C'(d) < oo 
such that if C > C , then for all x G W and R G [2(ji — 5), 2(p + 5)] we have 




ball 5(0, C)). Let 



f(C,R) 



C(B(0,C)nB(x,R)) 
C(B(x,R)) 




where R 1 = ^/R 2 - (1 + 2e) 2 and R 2 = ^/R 2 - (1 - 2e) 2 . 
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Proof. For z = (zi, z 2 , . . . z d ) G M. d we write z = (zx,z 2 ) and z = 
(z 3 , Z4, . . . z d ), and we write p and p for Euclidean distances on M 2 and M d ™ 2 
respectively. 

We claim that for any x G W, 

B 2 (w',e) x [B d _ 2 (0,C)nB d _ 2 (x,R 1 )] 

B 2 (w', e) x [B d _ 2 (0, C) n B d - 2 (x, i? 2 )] . 
To prove this, first note that for x G B 2 (w, e) and z G B 2 (w', e), we have 
it/) — p(z, w') — p(x, w) < p(z,x) < p(z, it/) + p(w, w') + p(w', x), 
which gives 

1 - 2e < p(z,x) < 1 + 2e. 

Take any 2 G W PI B(x, R). Then 2 G B 2 (w', e) and x G B 2 (w, e). Since 
z G B(x,R), we have 



pfoz) < \ R 2 ~ pU,x) 2 < ^R 2 - (1 - 2e) 2 = i? 2 



that is z G B d _ 2 (x, R 2 ). The second inclusion of the claim follows. 

On the other hand, if z G B 2 (w',e) x (B d _ 2 (0,C) H B d _ 2 (x,Ri)), then 
p(z,x) < l + 2e and p(z, x) < This gives p(z, x) < R, so 2 G W'nB(x, R). 
Thus the first inclusion of the claim is proved. 

The required statement now follows from the claim. For the lower bound 
we use Lemma |7J noting that 

R x G [ v / 4(/i-5) 2 -(l + 2e) 2 , v / 4(/i + 5) 2 -(l + 2e) 2 ] C (0,oo). 

For the upper bound we discard the intersection with B d _ 2 (0, C). □ 

Proof of Lemma O We have 

CS = C(W'nB(x,r + fi + 5)) - C(W' n B(x,r + p - 5)). 

Therefore, choosing C according to Lemma |H1 for C > C we have by that 
lemma, 



CS > us 2 u) d - 2 e 2 



r d ~ 2 /3 - r d ~ 2 
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where r := y/(r + /i + S) 2 - (1 + 2e) 2 and r := ^/(r + fi-5) 2 - (1 - 2e) 2 . 
For C chosen as above, we want to find a uniform lower bound on £S for all 
possible r £ [/i — 5, /i + 5]. We claim that for d > 11, the function 

^(r) :=r d - 2 /3-r Q! - 2 

is increasing on r £ |ju — 5, // + 5] . Once this is proved we obtain 

g(r) > gQi -S)> (Vh2) d - 2 /3 - (V^73) u ~ 2 > | (l.2^ - 1 

on this interval, and the result follows. 
To prove the above claim note that 

d-2 



g'(r) 



(r + /j, + <5)r d " 4 - 3(r + fi - 8)r 



d-4 



so it is enough to check that 



~2 I > 3 



for the required values of r and d. But is is straightforward to check that 
the quadratic function r 2 — 1.4r 2 is positive for r £ \p — 5, \i + 8], hence 
r 2 /r 2 > 1.4; and we have 1A^~ > 3. □ 



Proof of Theorem 2 



Proof of Theorem [21 Let d > 45, and construct the hard sphere process 
T as in Section |2J where the constant C is chosen according to Lemma and 
the intensity A of the Poisson process will be chosen later. 

Consider the first layer of the construction, and recall that at step n, 
vertex w n is explored, and a sphere of radius r Wn is placed with centre x Wn = 
Vw„ ^ W(w n ), provided the vertex is found to be good. As a notational 
convenience, if the algorithm stops during step N we write w n = A for all 
n > N, and call A a good vertex. Let T n be the a-algebra generated by all 
of the random variables 

( Wi, y Wi , X Wi , r Wi , n|vv(iui)5 n\B(x Wi ,r Wi ) ) ._ i 
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that is "the information known up to and including step n" . Let T-\ be the 
trivial a-algebra. 

We will compare the set of good vertices with a percolation cluster. Sup- 
pose that for some q we can show that for all n > 0, 

P(w n is good | Fn-i) > 1 almost surely (1) 

(that is, each newly explored vertex is good with probability at least q uni- 
formly in the past). Then the random set of good site vertices stochastically 
dominates C^(0), the open cluster at the origin for site percolation with 
parameter q 2 on H. This is because, as long as it is possible to add a new 
site vertex to the cluster of good site vertices at the origin, the algorithm 
attempts to do so by first exploring the intervening bond vertex, and then 
immediately exploring the new site vertex. The probability that both steps 
succeed is at least q 2 . Therefore we can compare with a cluster-growing 
algorithm for C*|(0). 

Also recall that the construction gives, in each layer, a cluster of touching 
spheres (with radii £ \p — 8, \i + 8]) including a sphere with its centre in 
the cell of each good vertex. Note therefore that if J is such a cluster of 
spheres then J^ c+< ^ contains all good vertices. Therefore, since the layers of 
the construction are independent, if we can establish with q 2 > p s c ltc (M) 
then the statements of the theorem will follow by Lemma El Thus all that 
remains is to prove (0) for with q > 0.892 > V0.794 > ^pf tc (E). 

Recall that in step n of the construction we choose a random Poisson 
point y = y Wn (if any exists) in a certain J-^i-measurable set S, and check 
to see whether a certain ball B(y, r Wn ) is free of other Poisson points. If both 
steps succeed then the vertex is declared good. The radius r Wn depends on 
which point y is chosen, but since it can be at most fi + 5 we can bound 
the required probability by the probability that the larger ball B(y,fi + 8) 
contains no other points: 

P(u>„ is good | J- n -i) > P(y ^ A, and y is (fi + 5)-isolated | F n -\). 

The latter event depends on II only through its restriction to W(w n )^ +S \ 
The conditioning on T n -\ does not affect the process n|^( Wn ), while the fact 
that one vertex v n adjacent to w n is good tells us only that a certain T n -\- 
measurable set Z not intersecting S contains no points of II. (Recall that 
2(n + 8 + e) < y/3, so the conditioning on non-adjacent vertices has no effect). 
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Thus we deduce by Lemmas and |U that 



P{w n is good | ^ n _i) > e ~ XCB{ ^ +5) - e~ xcs 

> 1 - XCB(0,fi + 6) - e~ xcs . 

But Lemma El applies to give 

almost surely. 

Thus, we require that 



F(X) := 1 - XB - e~ XA > 0.892 

where 



d-2 



A := 3 [1.2— - 1J and B := ^0.85". 

For each d we can choose the intensity A so as to get the best bound. Differ- 
entiating shows that F has a maximum at 

A* = A*(d) = ^log|, 



at which 



A/5 



But co> 2 = 7T and = ^ (see for example P), therefore 



A 10" 4 rf (1.2V - 1 
B ~ 6 x 0.85 d ' 

Thus A/B is an increasing function of d, and it is easy to check that for d > 31 
we have A/.B > 1, and therefore A* is positive (which is a requirement for an 
intensity). Furthermore, F(X*) is an increasing function of A/B for A/B > 1, 
and therefore an increasing function of d > 31. Finally it is straightforward 
to check that F(X*(d)) > 0.892 for d > 45, as required. □ 
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Remark — choice of lattice. Our construction could be adapted to work 
for other two-dimensional lattices such as the square lattice in place of the 
hexagonal lattice, but at the expense of increasing the dimension d. The 
fundamental requirement is that the diagonals (that is, the distances between 
non-adjacent vertices of the lattice) must be greater than 2/^/3 times the edge 
length. This allows the set S (in which we try to find a possible centre for 
a sphere) to be made much larger in volume than the hard spheres (which 
must be empty of other Poisson points), while preventing overlap between 
non-adjacent spheres. 

4 Stationarization 

In this section we deduce Theorem ^ from Theorem EJ 

Proof of Theorem ^ For each positive integer n, let U n be a random 
variable uniformly distributed on B(Q,n), and independent of T. Define 
the randomly shifted process T n = F + U n . Clearly Y n is a percolating 
Poisson hard sphere process, and has no spheres larger than K. We shall 
use Prohorov's Theorem to construct A as a weak limit of the sequence (T n ), 
and show that it has all the required properties. 

We claim that the sequence of random variables (T n ) is tight in the weak 
topology induced by the vague topology on point measures on M. d x [0, oo). 
To check this, it is enough to show that the sequence (T n (A)) is tight for any 
relatively compact Borel iGM^x [0, oo). (See [3] Lemma 16.15). Any such 
A is a subset of B x [0, oo) for some bounded Borel B C R d , and we have a.s. 
T n (^4) < r«(-B x [0, K]) since T n has no spheres larger than K. But the latter 
quantity has the same (Poisson) distribution for each n, so the sequence is 
clearly tight as required. 

Now let A be any weak subsequential limit of (r n ), so 

T nfe =>- A as k -> oo (2) 

in the topology referred to above. Clearly A is integer-valued and thus a 
sphere process. Furthermore, it is easily seen that the set of hard sphere 
processes supported a.s. on WL d x [0, K] is weak closed, and therefore A is a 
hard sphere process supported on M. d x [0,-fT]. Let B C M. d be a bounded 
Borel set with £-null boundary. By Theorem 16.16 of 0, the convergence in 
(j2J) implies the convergence in distribution T Hk (B x [0, K]) — > A(B x [0, K\). 
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Hence the latter has Poisson distribution with mean LB and so A is a Poisson 
sphere process. 

Next we show that A is invariant. It is sufficient to show that for any 



z g 



and any continuous compactly-supported function / : 



d x [0, oo) -> 
T + U n ; we 



[0, 1] we have Ej f dA = Ej f d(A + z). Recall that T n 
shall compare T n and T n + z for large n. Fix z and /, and write J = J{n) = 
-8(0, n)nB(z, n). Let V be random variable uniformly distributed on J. Also 
write a = a(n) = £(B(0,n) \ J)/£B(0,n), and note that a < cn d - x jn d -> 
as n — >• oo. Recall that U n was uniform on -8(0, n). Observe that conditional 
on U n G J (which occurs with probability 1 — a), the law of {7 n equals that 
of V. And conditional on U n + z G J (which occurs with probability 1 — a), 
the law of U n + z equals that of V. Hence 



E fdT n 



E 



U n E J 



1-aO+E 



U n e.P 



a 



= E J f d{T + V){l-a)+E Jfd{T + U n ) 
and similarly 

E [ f d(T n + z) 



U n eJ c 



a . 



E / f d(T + V){l-a) + E 



f d(T + U n + z) 



U n + z G J c 



a. 



Since U n is independent of T, the processes Y + U n and r + U n + z appearing 
on the right sides of the above equations are Poisson sphere processes, even 
when conditioned on U n . Now for any Poisson sphere process T say, E J fdT 
is bounded by the expected number of Poisson points in the projection of 
the support of / onto R d , that is the Lebesgue measure of that projection, 
C say. Hence, subtracting the two equations above gives 



E / / dT n — Elf d{T n + z 



< + aC — > as n — ► oo. 



Taking weak limits as k — >• oo of T nk , T nk + z we deduce from (J2J) that 



E / / dA - E / / d(A + z) 



17 



Hence A is invariant as required. 

Finally we must show that A percolates almost surely. For a hard sphere 
process T and for < a < b let H a ^(T) be the event that G(T) has a 
connected set of spheres with radii at most K which intersects both -6(0, a) 
and B(0,b) c . Also let H at00 (T) be the event that G(T) has an unbounded 
connected set of spheres with radii at most K which intersects -8(0, a). Note 
that H a oo {T) is the decreasing limit of the events H ab {T) as b —> oo. Also 
denote by I(T) the union of all infinite clusters of the hard sphere process 
T. Recalling the definition of T n above, we have 

P(H a>b (v n )) > p(# ai00 (r n )) 
= P(oe/(r)W) 

£(/(r)W n g(0,n)) 
£B(0,n) 

Hence by Fatou's Lemma and the definition of lower density, 
liminfP(# a)6 (r n )) > EP(/(r)W). 

Note that the event H a ^(T) depends only on the process T restricted to the 
compact set B(0,b + K) x [0, iT]. Furthermore it is straightforward to see 
that the event is closed in the vague topology on point measures. It follows 
from (J2J and the Portmanteau Theorem ([3] Theorem 4.25) that 

P(# a , 6 (A)) > limsu P P(if a>6 (r re J). 

From the last two inequalities we have 

P{H a>b (A)) > E2?(/(r)W), 
and hence letting b — > oo we deduce 

P(A percolates) > P(H a>00 (A)) > EV(I(T) {a} ). 
Finally letting a — > oo, Theorem |21 gives 

P(A percolates) = 1. 

□ 
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Remark — positive radii. As noted earlier, our proof may be adapted 
so that all the spheres have positive radii. To achieve this, take a small 
parameter rj > 0, and modify the construction in Section |21 as follows: having 
chosen a potential centre y w for a sphere, we declare the vertex good only if 
the larger ball B(y w ,r w + 77) contains no other Poisson points (rather than 
the ball B(y w ,r w )). If 77 is small enough then this does not affect the later 
computations, and we still obtain a non-invariant percolating hard sphere 
process V for d > 45. But now T has the additional property that no zero- 
radius sphere is within distance 77 of any nonzero-radius sphere. Therefore 
the stationarized version A will inherit the same property. Finally, we modify 
A as follows. If there is a zero-radius sphere centred at z, replace it with a 
sphere of radius r, where r is 1/2 of the distance from z to the closest other 
sphere of A (including other zero-radius spheres). This r is always positive 
because the radii of the existing spheres are uniformly bounded above, and 
any bounded region of lR d contains only finitely many Poisson points. 

Open Problems 

(i) Does there exist a percolating Poisson hard sphere process (invariant 
or non-invariant) in dimensions 2 < d < 44? The case d = 2 seems 
particularly interesting. 

(ii) In any dimension, does there exist a percolating, invariant Poisson hard 
sphere process which is a deterministic function of the Poisson process? 

(iii) Do percolating hard sphere processes exist for other point processes, 
such as Gaussian zeros processes 
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